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Abstract — Trajectory  generation  for  quadruped  robots  is  a 
challenging  task  since  they  are  underactuated  systems  which 
must  balance  using  sensory  feedback  and  satisfy  ground  contact 
constraints.  There  is  a  substantial  body  of  evidence  that  many 
animals  use  central  pattern  generators  (CPGs)  for  generating 
joint  trajectories  and  regulation  through  sensory  feedback. 
However,  CPG  models  formulated  in  the  joint  space  do  not 
explicitly  formulate  or  account  for  ground  contact  constraints, 
especially  during  turning  which  can  introduce  foot  slip  in  the 
gait.  Task-based  CPGs  offer  several  advantages  in  that  they  can 
explicitly  satisfy  ground  contact  constraints,  and  we  suggest  an 
approach  to  generating  foot  velocity  controls  in  the  body  frame 
to  enable  omnidirectional  locomotion. 

I.  Introduction 

Legged  robots  such  as  quadrupeds  have  great  potential  for 
navigating  rough  terrain  that  other  autonomous  vehicles  can¬ 
not.  Many  animals  are  able  to  quickly  negotiate  rough  terrain 
such  as  steep  mountain  slopes.  However,  legged  robots  are 
often  relegated  to  the  confines  of  a  laboratory  setting  where 
either  the  ground  is  flat,  solid,  and  devoid  of  obstacles,  or 
the  environment  is  known  beforehand  and  accurate  global 
pose  information  is  available.  Motion  planning  techniques 
can  be  used  to  generate  collision  free  paths  when  the  envi¬ 
ronment  is  known  in  advance.  To  reduce  the  dimensionality 
of  the  planning  problem,  a  high  level  motion  planner  can 
generate  paths  for  the  base  of  the  robot  to  follow,  while  foot 
placements  are  planned  online  which  maximize  the  stability 
margin  from  the  zero  moment  point.  This  approach  was  used 
effectively  on  the  Little  Dog  platform  [1],  However,  such 
approaches  assume  the  robot  has  access  to  terrain  maps, 
which  are  constructed  offline  and  require  access  to  global 
pose  information  by  means  of  a  motion  capture  system. 

To  get  these  robots  out  of  the  laboratory  and  into  unknown 
field  environments,  planning  collision  free  paths  without 
external  sensing  is  of  utmost  importance.  Large  uncertainties 
in  sensing  capabilities  and  especially  global  pose  make  it 
difficult  to  perfectly  negotiate  obstacles  in  a  collision  free 
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Fig.  1 .  The  Allegro  Dog  quadruped  platform 


manner.  As  Kalakrishnan  et  al.  note,  compliance  is  critical 
for  interaction  in  uncertain  environments  [1],  Further,  reflex 
behaviors  are  effective  in  providing  quick  stabilizing  maneu¬ 
vers  when  interacting  with  unknown  terrain  [2].  In  order  to 
navigate  in  rough  terrain  in  GPS-denied  environments,  we 
propose  having  the  robot  base  follow  a  desired  collision  free 
path  which  is  planned  in  local  coordinates  while  footholds 
are  generated  reactively  using  a  central  pattern  generator 
(CPG).  This  approach  requires  omnidirectional  locomotion 
(the  ability  to  start,  stop,  and  move  in  any  direction  with 
a  desired  twist  imparted  on  the  body)  should  the  vehicle 
deviate  off  course  and  must  return  to  its  intended  path. 

There  is  substantial  evidence  that  many  animals  (mam¬ 
malian  quadrupeds,  bipeds,  lamprey)  have  CPGs  in  the  spinal 
cord  that  assist  in  the  generation  of  gaits  [3],  CPGs  have 
been  used  frequently  in  robotics  [4],  for  quadruped  robots 
in  particular,  to  generate  joint  trajectories  (gait  sequences) 
and  incorporate  sensory  feedback  to  stabilize  locomotion  [5], 
[6]  .  However  (with  notable  exceptions  [7],  [8]),  the  bulk 
of  research  has  focused  on  design  of  steady  state  forward 
gaits.  The  majority  of  CPG  models,  when  they  do  specify  a 
turning  gait,  execute  turning  maneuvers  by  changing  the  joint 
bias  of  the  hip  joints  requiring  specialized  arm  geometry. 
This  approach  has  the  further  limitation  that  it  does  not 
explicitly  formulate  ground  contact  constraints,  which  can 
introduce  foot  slip  in  the  gait.  While  most  CPGs  generate 
joint  trajectories  directly,  Barasuol  et  al.  have  demonstrated 
using  task  space  CPG-like  coupled  oscillators  to  control  the 
HyQ  platform  [9].  The  result  is  that  fewer  oscillators  are 
required  and  direct  approaches  exist  for  enforcing  ground 
contact  constraints. 
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II.  Hardware  Description 


This  work  used  the  Allegro  Dog  robot  which  is  a 
commercial-off-the-shelf  quadruped  robot  from  SimLabs 
Co.,  Ltd.  (see  Figure  1).  The  system  is  450mm  in  length, 
240mm  wide,  stands  at  500mm  in  height  and  weighs  20kg. 
Each  leg  has  actuated  revolute  joints  for  the  hip-roll,  hip- 
pitch,  and  knee-pitch  degrees-of-freedom  and  a  passive 
spring-loaded  prismatic  joint  below  the  knee  (totaling  to 
twelve  actuated  and  four  passive  joints).  Each  revolute  joint 
is  controlled  with  a  90W  Maxon  brushless  DC  motor  with  an 
86:1  ceramic  planetary  gearhead  and  is  sensed  directly  at  the 
motor  with  a  1000  counts  per  turn  (CPT)  incremental  mag¬ 
netic  encoder.  Homing  switches  are  used  during  initialization 
to  find  the  mechanical  zero  position  of  each  joint.  Motor 
controllers  close  a  current  loop  at  1kHz.  Encoder  values  are 
sent  to  and  current  commands  are  received  from  the  Allegro 
Dog  embedded  controller  over  a  CAN  bus  at  500  Hz. 

The  base  body  has  a  CH-Robotics  UM6  inertial  mea¬ 
surement  unit  (IMU)  that  provides  body  orientation,  angular 
velocities,  and  linear  accelerations.  This  IMU  updates  inter¬ 
nal  state  at  500  Hz  and  transmits  data  to  the  Allegro  Dog 
embedded  controller  at  about  70  Hz  through  a  TTL-USB 
converter. 

Power  is  provided  externally  through  a  tether,  or  onboard 
with  24V/4500mAh  batteries.  The  Allegro  Dog  embedded 
controller  has  both  wired  and  wireless  Ethernet  capabilities. 
All  commands  to  and  state  updates  from  the  Allegro  Dog 
embedded  computer  are  communicated  over  Ethernet  using 
the  Lightweight  Communications  and  Marshalling  LCM 
library  [10]. 

III.  Task-based  Central  Pattern  Generator 

We  use  a  task-based  CPG  based  which  is  an  extension 
of  the  coupled  oscillators  in  [9].  This  CPG  consists  of 
a  network  of  coupled  modified  Hopf  Oscillators  and  a 
filtered  system  output  that  can  be  used  to  directly  modify 
the  foot  ground  velocities  when  in  stance  mode.  Define 
Xj  =  ( Xi,yi,Zi )  G  R3  to  be  the  task  position  of  the  ith 
foot  where  the  total  number  of  feet  is  N  =  4.  Combining 
all  the  task  positions  into  a  single  3  x  N  matrix  yields 
X  =  [xi,---  , xjv]  Let  xpo  i  be  the  center  position  of  the 
ith  foot  oscillator’s  limit  cycle.  The  displacement  about  the 
center  position  is  denoted  as  x,  =  x,  —  xpo  ,:  +  Axt.  In  this 
case  AX  =  [Axi,  •  •  •  ,  Axjv]  =  0.  The  equations  of  motion 
of  the  CPG  are  in  (1): 
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Define  the  instantaneous  phase  dependent  angular  fre¬ 
quency  iOi  depending  on  the  leg  swing/stance  phase.  Also 
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Fig.  2.  Limit  cycle  of  an  example  foot  position  controlled  by  the  CPG. 
The  solid  black  line  is  the  limit  cycle  of  the  canonical  CPG  (1)  while 
the  gray  dotted  line  the  filtered  output  (2),  which  is  smoothly  clipped  at 
z  =  Ztd,i  =  —0.1  Hs. 

note  parameters:  gains  ct,/3,7  >  0,  stride  length  Ls,  step 
height  Hs,  swing/stance  duty  factor  Df  and  desired  forward 
velocity  Vfjes.  Cy  is  the  coupling  matrix  defining  foot 
phase  relationships  and  are  defined  in  [9],  Gaits  such  as 
walk,  trot,  and  bound  can  be  characterized  by  Df  and 
C ij.  A:  is  a  scalar  parameter  to  regulate  the  strength  of  the 
coupling.  The  gating  functions  crcl  i(zc  i)  =  ( e~bZi  +  1)_1 
and  (Jc‘2.i(zc,i)  =  ( ebZi  +  1)_1  are  sigmoids  which  enable  a 
smooth  transition  between  the  swing  angular  frequency  and 
stance  angular  frequency.  Note  that  the  z  component  is  used 
to  identify  whether  the  foot  is  in  swing  or  stance  mode. 

The  filtered  output  dynamics  are  as  follows: 

Nf (xC)i  T  iTc(xc  ?;  Xy  .j))  (Tyi  i  (zc 

-ViO  f2,i{zc,i)  (2) 

vn,i(zc,i)  =  {e-b^'i-Ztd’i)  +  I)”1 
=  (e6^-****)  +  I)”1 

The  use  of  the  filtered  output  and  gating  functions  <7/1  , 
and  <Jf2,i  permit  the  smooth  mixing  of  different  task  kine¬ 
matics  for  whether  the  leg  is  in  swing  or  in  stance  mode. 
When  in  swing  mode,  the  filter  output  tracks  the  canonical 
CPG.  When  the  feet  are  in  stance  mode,  the  task  kinematics 
are  defined  by  V)  =  [VitX,Vi>y,VitX]  G  R3  which  can  be 
computed  to  enable  omnidirectional  locomotion.  The  step 
depth  parameter  zt(u  is  a  controllable  input  for  smoothly 
clipping  the  ellipse  (See  Fig.  2).  However,  one  must  be 
careful  to  ensure  that  the  canonical  CPG  (1)  and  filtered 
output  (2)  trajectories  are  well  matched  to  minimize  aberrant 
behavior. 

IV.  Extensions  to  the  CPG  to  enable 
Omnidirectional  Locomotion 

We  extend  the  CPG  defined  in  (1)  and  (2)  to  be  capable 
of  omnidirectional  locomotion.  We  modify  the  dynamics  of 
the  CPG  to  generate  smooth  trajectories  that  execute  the 
following  actions: 

1)  The  ability  to  start  and  stop.  This  requires  adding 
an  additional  control  input  which  adds  a  supercritical 
bifurcation  point  in  the  Hopf  Oscillator.  See  Sec.  IV-A. 


2)  Change  the  direction  of  the  angular  velocity  to  enable 
forward/reverse  locomotion.  See  Sec.  IV-B. 

3)  The  ability  to  march  in  place  for  load  testing,  turn  in 
place  maneuvers  and  diagnostic  purposes.  See  Sec.  IV- 
B. 

4)  Automatically  generate  turning  and  crab  gaits  by  gen¬ 
erating  a  body  twist  controller  to  generate  foot  veloc¬ 
ities  that  reproduce  the  desired  locomotion.  See  Sec. 
IV-C. 

The  modified  dynamics  are: 
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Where  p,r]  £  {—1,1}  are  control  inputs  to  regulate 
locomotion  behavior.  Filter  dynamics  are  only  modified  for 
march-in-place. 


A.  Start  and  Stop 

Previous  work  [7]  has  modified  the  Hopf  oscillator  to 
contain  a  supercritical  Hopf  bifurcation,  which  enables  the 
oscillator  to  enter  two  different  modes: 

•  p  =  —  1:  discrete  motion  to  a  user-selectable  stable 
equilibrium  point 

•  p  =  1:  rhythmic  locomotion  tending  to  a  limit  cycle 

The  effect  of  switching  p  on  the  oscillator’s  performance 

can  be  seen  in  Fig.  3.  However,  merely  modifying  the 
canonical  CPG’s  dynamics  is  not  sufficient  to  have  start/stop 
gait  transitions  as  tracking  issues  will  occur  between  the 
canonical  system  and  the  filtered  output. 

In  addition,  AX  can  be  used  to  change  the  stable  equilib¬ 
rium  point  of  the  oscillator.  When  in  stop  mode,  AX  should 
be  selected  such  that  Xpo  +  AX  is  within  the  attractive  basin 
and  will  quickly  converge  to  the  limit  cycle  of  the  desired 
gait.  Otherwise,  AX  =  0. 


oscillator  changes.  One  can  imagine  this  being  the  result 
of  shrinking  the  diameter  of  a  wheel  in  a  differential  drive 
robot.  Despite  maintaining  angular  rate,  the  ground  velocity 
will  decrease.  By  observing  that  along  the  limit  cycle  of  the 

oscillator  when  p  =  1,  £%’*-! =  1,  one  can  scale  Vftdes 

and  estimate  the  CPG’s  instantaneous  forward  velocity  Vf  of 
the  canonical  CPG  given  the  current  state,  then  average  the 
result  of  all  the  legs  so  that  the  canonical  CPG  and  filtered 
output  are  matched.  See  (4): 


Fig.  4  demonstrates  how  setting  Vi  =  [—  Vf  dest  0, 0] 
causes  tracking  issues  between  the  canonical  CPG  and  the 
filtered  output,  while  Vi  =  [—  V/,0,  0]  improves  tracking 
performance  of  the  filtered  output. 


(a)  Having  fixed  Vf  during  startup  and  shutdown  sequences  causing 
significant  tracking  issues  in  the  filtered  output. 
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(b)  Proper  scaling  of  Vf  using  (4)  enables  the  filtered  output  to 
properly  track  the  canonical  CPG 

Fig.  4.  Tracking  issues  between  the  canonical  CPG  and  filtered  output 
arise  during  startup  and  shutdown  sequences.  Using  (4)  to  account  for  the 
instantaneous  forward  speed  of  the  CPG  corrects  for  this  effect. 


Fig.  3.  The  x  component  of  a  Hopf  oscillator  changing  between  start  and 
stop  modes  by  switching  control  input 

While  angular  velocity  is  preserved  in  the  Hopf  oscillator 
during  shutdown,  the  translational  velocity  with  respect  to 
the  ground  plane  is  not  preserved  as  the  amplitude  of  the 


B.  Forward/Reverse  and  March  In  Place 

To  switch  between  forward  and  reverse  locomotion  by 
modifying  the  original  canonical  CPG  (1),  define  angular 
direction  r/  such  that  77  =  1  for  forward  locomotion  and 
77  =  — 1  for  reverse  locomotion. 


World  frame  arbitrary 


Fig.  5.  Schematic  representation  of  the  coordinate  frames  and  the  virtual 
rigid  body  used  in  the  foot  velocity  controller.  Imparting  a  twist  or  spatial 
velocity  on  the  robot  base  will  induce  motion  of  the  feet  and  must  be  actively 
canceled  by  the  foot  controller  to  prevent  slip. 


It  is  also  possible  to  redefine  the  filter  dynamics  to  enable 
a  simple  march  in  place.  One  can  replace  the  x component 
of  the  filter  dynamics  to  have  stable  linear  dynamics  decaying 
to  the  foot  center  position  Xpo,i  when  the  leg  is  in  swing 
mode.  Summarizing  (2)  as  x =  [/i,  /2,  /j]\  augment  filter 
dynamics  to  be  (5): 
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Where  d  £  V  is  a  new  control  input  specifying  the 
hybrid  state  of  the  system  (walk,  stop,  march,  etc),  which 
is  explained  further  in  Sec.  IV-D.  Passing  V i  through  the 
modified  march  dynamics  (5)  when  the  leg  is  in  stance  mode 
can  allow  the  quadruped  to  rotate  in  place  when  using  body 
twist  control  (See  Sec.  IV-C). 


C.  Foot  Velocity  Generation  for  Body  Twist  Control 

Suppose  the  base  of  the  robot  is  expected  to  generate 
motion  constrained  to  a  plane  parallel  to  the  ground  plane 
with  a  given  arbitrary  twist  command  |  £  se( 2),  with  twist 
coordinates  £  =  (£x>  €y,  &,),  where  is  the  desired  forward 
base  velocity,  £y  is  the  side  step  velocity  (for  crab  walk), 
and  is  the  desired  turning  rate  along  the  z  axis  (upward). 
Using  these  coordinates,  the  twist  is  in  (6): 
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Note  that  se(2)  is  the  Lie  algebra  of  SE(2)  so  is  a 
homogeneous  transform  £  SE( 2)  [1 1].  It  is  assumed  that  this 
twist  command  is  given  in  the  robot  body’s  frame  instead  of 
any  particular  spatial  frame  (or  inertial  frame).  The  position 
of  the  foot  should  not  move  in  the  world  frame  or  any 
spatial  frame  so  the  footholds  do  not  slip  by  controlling 


the  position/velocity  of  the  feet  with  respect  to  the  body. 
One  spatial  frame  of  interest  is  a  spatial  frame  positioned 
at  the  foot  and  oriented  as  the  robot  base  frame  in  Fig. 
5.  However,  consider  a  virtual  rigid  body  where  the  point 
contacts  of  the  feet  in  stance  mode  are  in  a  fixed  position  with 
respect  to  the  body  frame  b.  The  movement  of  the  base  will 
induce  motion  of  the  point  contacts  of  the  stance  feet  since 
they  are  a  point  on  the  rigid  body.  Further,  assume  the  current 
task  positions  of  the  feet  are  fully  inside  the  workspace  of 
the  robot  and  motions  in  all  directions  are  admissible. 

Define  b  £  se(2)  to  be  the  spatial  velocity’  of  the  base 
moving  with  respect  to  the  spatial  frame  placed  at  foot  i 
which  is  in  stance  mode  (7).  The  desired  body  velocity  V*b 
can  be  transformed  into  the  spatial  velocity  V£b  using  the 
adjoint  transform  [11]: 
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6- 
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0  £y  + 
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(7) 


For  there  to  be  no  slipping  of  the  stance  feet,  we  must 
counter  the  spatial  velocity  of  the  feet  using  control  input  Vt 
in  (2)  such  that  the  no  slip  constraint  (8)  holds: 


([ii  m  o] V  =  v:ibp°  -  Vi  =  0  (8) 

where  p?  =  [.T;,  y{,  if  =  [0,0,1]'  is  the  origin  of  spatial 
frame  Wi  placed  at  the  foot  location  and  v,  =  [«,>,  Vi.y,0]\ 
The  quantity  V^.^pf  is  the  linear  velocity  of  the  point  pf 
induced  by  spatial  velocity  V£  b  in  the  spatial  frame.  Note 
that  v,  and  (8)  are  in  homogeneous  coordinates  while  Vt  £ 
M3  where  we  have  the  planar  assumption  ij  =  Vl  z  =  0. 
Satisfying  (8)  requires  that  CPG  control  input  V,  be  of  the 
following  form  (9): 


r^i,x 

£x  Vf,i£uj 

= 

vi,y 

= 

0 

0 

Figures  6  and  7  highlight  the  task  trajectories  of  the 
feet  in  the  body  frame  when  executing  a  turning  maneuver. 
Considering  the  previously  observed  tracking  issues  in  the 
start/stop  motion  primitives  in  Sec.  IV-A,  the  desired  twists 
are  normalized  with  respect  to  the  instantaneous  forward 
velocity  ||(£x,£!/)||  =  Vf  so  that  both  the  canonical  CPG  and 
filtered  output  are  closely  matched  during  startup/shutdown 
sequences.  This  technique  will  generalize  to  twist  controls 
in  se(3). 

D.  Summary:  Hybrid  State  Definitions 

We  found  that  in  modifying  a  CPG  based  gait  controller, 
it  is  helpful  to  have  the  following  hybrid  state  d  £  V 
to  control  the  rhythmic  aspects  of  the  CPG  where  V  - 
{dgtop i  d f  oriuardi  dreversei  dmarch} •  This  enables  the  CPG 
to  generate  smooth  transitions  to  start  and  stop  walking 
without  explicitly  encoding  transition  sequences.  Parameter 
configurations  for  each  hybrid  state  are  summarized  in  Table 
I.  Note  that  the  hybrid  state  is  not  assigned  a  particular  body 
twist  £.  Each  hybrid  state  can  reproduce  a  variety  of  stable 
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Fig.  6.  Example  reference  foot  trajectories  generated  by  the  CPG  when  a 
turn  command  is  issued  f  =  (Vy  des,0,ujjeyt).  The  desired  body  twist  is 
displayed  (linear  component  as  a  vector  and  rotational  component  as  an  arc) 
and  foot  task  trajectories  are  shown  in  the  body  frame  with  a  top-down  view. 
Note  that  the  Ackermann-like  steering  geometry  emerges  automatically  from 
explicitly  satisfying  ground  constraints. 


gaits  for  a  particular  set  of  £  €  se( 2).  This  set  is  dependent 
on  CPG  parameters  such  as  Xpo,  Ls  and  Hs.  For  a  quasistatic 
walk  gait,  criteria  such  as  zero  moment  point  could  be  used 
to  verify  the  feasibility  of  a  body  twist  for  a  given  hybrid 
state. 


state 

dstop  d  forward  d reverse  ^marc/i 

(v,v) 

AX 

{-l,Vprev.)  (1,1)  (1,-1)  (1,1) 

See  Sec.  IV-A  0  0  0 

TABLE  I 

A  SUMMARY  OF  HOW  THE  HYBRID  STATE  CONFIGURES  CPG  CONTROL 
VARIABLES. 


V.  Experimental  Setup  and  Results 

Since  the  Allegro  Dog  does  not  have  the  requisite  sensors 
to  be  able  to  measure  or  estimate  the  body’s  motion  with 
respect  to  the  ground  plane,  the  body  frame  of  the  Dog 
was  tracked  using  an  external  motion  capture  system  (Vicon) 
which  measures  the  Dog’s  pose  at  about  40  Hz.  This  data 
was  in  turn  synchronized  to  LCM  state  messages  regarding 
desired  and  actual  joint  trajectories,  discrete  motion  com¬ 
mands  (forward,  reverse,  etc)  and  the  control  twist  which  is 
internally  generated  by  the  foot  velocity  controller. 

While  the  Vicon  system  returns  poses  in  SE( 3),  it  would 
be  of  interest  to  extract  instantaneous  twists  between  poses. 
However,  to  numerically  differentiate  the  pose  directly  would 
be  highly  sensitive  to  the  noise  present  Recursive  Bayesian 
weighted  regression  is  used  for  automatic  (non-Gaussian) 
outlier  detection  with  a  forgetting  factor  A  =  0.95  [12], 

The  Allegro  Dog  uses  a  joint  space  trajectory  tracking 
PD  controller  with  fixed  base  (differential)  closed  loop 
inverse  kinematics  [13]  and  with  base  attitude  compensation 
(kinematic  adjustment  in  [9]). 

Figure  8  compares  the  control  twist  starting  from  standstill 
to  steady  state  forward  locomotion  £  =  (0.11®,  0,0),  and 
the  measured  twist  (angular  rate  from  IMU,  linear  velocity 
from  Vicon)  in  the  body  frame.  Unmodeled  dynamics  of  the 
foot-ground  contacts  cause  the  platform  to  sway,  contributing 
to  errors  in  the  reproduction  of  the  desired  twisL  Notice 
the  lateral  sway  in  vy  during  startup.  The  oscillations  in 
the  yaw  rate  uiyaw  are  unbiased  so  the  robot’s  net  motion 
over  time  is  forward  with  oscillations  in  yaw  angle  being 
less  than  3°  peak  to  peak  during  the  gait  cycle.  Further,  the 
forward  velocity  Vj^ea  =  0.11®  is  normalized  with  respect 
to  the  current  state  of  the  CPG  using  (4)  hence  the  minor 
oscillations  in  vx. 


VI.  Conclusions 


Fig.  7.  Example  reference  foot  trajectories  generated  by  the  CPG  when  a 
left  turn  command  is  issued  £  =  (Vf  des,  0,uieft).  This  is  viewed  from  a 
perspective  view. 


A  task-based  CPG  has  been  proposed  that  is  capable 
of  omnidirectional  locomotion,  including  start/stop  behavior 
and  a  foot  velocity  controller  that  generates  stance  feet 
velocities  to  impart  any  desired  twist  on  the  robot  base. 
Preliminary  results  show  that  the  foot  controller  is  able  to 
impart  the  desired  twist  on  the  body. 


twist  command 


_ i _ i _ i _ i 

94  96  98  100 
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Fig.  8.  Results.  Twist  command:  standstill  to  forward  twist  command 
£  =10. 1 1  ®,0,0))  (solid  thick  black),  and  measured  twist  (dotted  blue/gray). 
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